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Abstract 

The transverse mass spectra of D, J/ip and ip' in Pb+Pb collisions at 158 A-GeV are 
studied within a hydrodynamical model of the quark gluon plasma expansion and 
hadronization. The model reproduces the existing data with the common hadroniza¬ 
tion parameters: temperature T = Tn = 170 MeV and average collective transverse 
velocity vt — 0.2. 

The possibility of observing the quark gluon plasma (QGP) in nucleus-nucleus (A+A) 
collisions motivated an extensive experimental program. Rich experimental results on 
hadron spectra and multiplicities are now available for a broad range of collision energy 
and for various colliding systems [jlj . The anomalies observed j2[ in the energy dependence 
of pion and strangeness production indicate || that in central collisions of heavy nuclei 
the threshold for QGP creation during the early stage of the reaction is in the region of 
low CERN SPS energies 40 A-GeV). 

The equilibrium hadron gas (HG) model describes remarkably well the hadron mul¬ 
tiplicities measured in A+A collisions at top SPS [|] and RHIC || energies, where the 
creation of QGP is expected. The extracted hadronization temperature parameter is sim¬ 
ilar for both energies T# = 170 ± 10 MeV. This is close to an estimate of the temperature 
Tc for the QGP-HG transition obtained in Lattice QCD simulations at zero baryonic 
density (see e.g. ||). One may therefore argue that the QGP created in high energy 
heavy ion collisions hadronizes into an (approximately) locally equilibrated HG and the 
chemical composition of this HG is weakly affected by rescattering during the expansion 
of the hadronic system [[7|. 

In our previous paper || we formulated the hypothesis that the kinetic freeze-out of 
J/'ip and ip' mesons takes place directly at hadronization and that those mesons there¬ 
fore carry information on the flow velocity of strongly interacting matter just after the 
transition to the HG. Based on the measured J/ip and ip' spectra in Pb+Pb collisions 
at 158 A-GeV |J and using the hypothesis of the statistical production of charmonia 
at hadronization (E3- E5- HI we extracted a mean transverse collective flow velocity of 
hadronizing matter: vt — 0.2. 
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The effect of the rescattering in the hadronic phase was recently studied within a 
“hydro + cascade” approach [|13|, [L4|. A+A collisions are considered there to proceed in 
three stages: hydrodynamic QGP expansion (“hydro”), transition from QGP to HG and 
the stage of hadronic rescattering and resonance decays (“cascade”). The change from 
“hydro” to “cascade” modelling takes place at T = Tc, where the spectrum of hadrons 
leaving the surface of the QGP HG transition is taken as input for the subsequent cascade 
calculations. The results of Refs. 0, [hi suggest that the transverse momentum (p T ) 
spectrum of hi may be weakly affected during the cascade stage even for central Pb+Pb 
collisions at the top SPS energy. This is because of the small hadronic cross section and 
large mass of the hi hyperon |R|. The corresponding calculations for charmonia are not 
yet performed within this model. 

Thus we are faced with an intriguing problem: if the above considerations for charmo¬ 
nia |1 and hi JP| Q are correct, their p T spectra should be simultaneously reproduced 
using the same hydrodynamic parameters, T # and u T . In this letter we demonstrate that 
such a description is indeed possible. The transverse mass tut {itlt = \jvr + m2 > where 
m is a particle rest mass) spectra around midrapidity in Pb+Pb at the SPS (158 A-GeV) 
were recently measured for hi by WA97 Collaboration [16, 17| and for J/ip and -0' by 
NA50 Collaboration 0. These spectra will be the subject of the present analysis. 

Assuming kinetic freeze-out of matter at constant temperature T, the transverse mass 
spectrum of i-th hadron species (with mass mQ in cylindrically symmetric and longitudi¬ 
nally boost invariant fluid expansion can be approximated as |18 


dNi f R , m T cosh y T 
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where yx = tanh -1 Vx is the transverse fluid rapidity, R is the transverse system size, K\ 
and J 0 are the modified Bessel functions. The spectrum ([[]) is obtained under assumption 
that the freeze-out occurs at constant longitudinal proper time r = sjft — z 2 , where t is 
the time and z is the longitudinal coordinate. Thus the freeze-out time t is independent 
of the transverse coordinate r 


The analysis of the numerical calculations of Ref. 


shows that the latter is approximately fulfilled. The quality of the approximation made 
gets better for considered here heavy particles and small transverse flow velocities because 
a possible deviation from Eq. ([[]) is proportional to p^vx/{^rn T T) and thus it decreases 
with increasing particle mass at constant px- 

In order to calculate ([!]) the function yx{f) has to be given. A linear flow profile, 
yx(j) = yix ax • r/R, is often assumed in phenomenological fits [R|. The numerical cal¬ 
culations of Ref. 0 justify this assumption. For heavy hadrons analysed in this work 
the condition nrii » T is always satisfied and, therefore, the asymptotic form for large 
arguments of K\ (x) ~ x~ l d 2 exp(— x) can be used in Eq. ©• At SPS energies typical 
values of vt are small (uf. << 1 ) and consequently cosh y T = 1 + r.v'Jp and sinh y T = 


The experimental my-spectra are usually parametrised by a function])]: 
dNi 
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1 It corresponds formally to neglecting the transverse flow (vt = 0 in Eq. ( 0 )) , but introducing instead 

an “effective” temperature. 
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where the inverse slope T* is extracted from the fit to the data. However, when Eq. (2|) 
is considered as an approximation of Eq. flU), the inverse slopes T* should depend on 
both rrii and px- The limiting cases of T* behaviour at low and high px can be easily 
studied using the small and large argument asymptotic of the modified Bessel function 1$ 
in Eq. (Q): 
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where the average velocity Vx in Eq. (|3|) is defined as v\ = f 0 R rdrv^(r )/ f R rdr. The 
maximum velocity v™ ax in Eq. (||) is related to vx as = (vx ax ) 2 /% provided a linear 
flow profile is assumed, Vx(r) = v™ ax ■ r/R. Note that T* in Eq. (3|) is equivalent to the 
well known “blue shifted” temperature, T[(l + v™ ax )/(l — r;™ 01 )] 1 / 2 (see e.g. |L4], |I8|), 
calculated for (v™ ax ) 2 << 1. The shape of the “high-p T ” tail (p T >> rrii) of the m T 
distribution © is “universal”, i.e. T* given by Eq. ([|) is independent of particle mass 
rrii. On the other hand, the inverse slopes T* ([3j) at “low -px" are strongly dependent on 
rrii. Two remarks are appropriate here. First, for heavy particles like H and J/ip the term 
\mpv\/T in Eq. ([|) is not small compared to one thus the second (approximate) equality 
in this equation is violated. Second, a condition of the validity of Eq. (||), pxv t « T, is 
too restrictive for heavy hadrons, e.g. for T = 170 MeV and vt — 0.2 discussed below it 
leads to m T — rrii « 0.3 /rrii GeV 2 /c 4 . This means that Eq. (j3|) is valid for the values of 
rrix — rrii which are much smaller than 0.2 GeV/c 2 for and than 0.1 GeV/c 2 for J/^. 

Summarising: (a) none of the asymptotic regimes (|3j) or (|]) can be clearly seen in the 
experimental mx spectra, i.e. neither “low-pr” (my — rrii « 0.3 /rrii GeV 2 /c 4 ) nor “high- 
Pr” (w-t — rrii » rrii) approximations are useful ones (at least for studying the available 
rrir spectra of and charmonia); (b) fitting the experimental mx spectrum of i-tli hadron 
species by Eq. (D) one finds in fact the “average inverse slopes” which depend not only 
on particle mass m,, but also on the rrix — rrii interval covered in a given experiment (see 
also Ref. Jl4|, where T* have been discussed separately for mx — m* < 0.6 GeV/c 2 and 
for 0.6 GeV/c 2 < m T — m* < 1.6 GeV/c 2 ). 

For small values of Vx relevant for our discussion a good approximation of Eq. ([[]) 
at mx — mi < m, can be obtained by substituting the vx distribution in Eq. (1|) by its 
average value Vx and by using large argument K 1 asymptotic: 
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We checked numerically that the values of the parameter Vx extracted from the fits (see 
below) to Eqs. ([!]) and (|5|) assuming a linear velocity profile differ by about 5%. 

We turn now to the test of our hypothesis of the kinetic freeze-out of J/4>i '0 / and 
occurring directly at hadronization i.e. at T = Th = 170 MeV. The m^-spectra 
of these hadrons are measured around midrapidity [[|. [l(| for Pb+Pb collisions at 158 
A-GeV. The fit to these data performed using Eq. (151) with T = T H = 170 MeV yields 
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vt = 0.194 ± 0.017 and x 2 /dof = 1.3. The value of vt varies by =f 0.016 when Th 
is changed within its uncertainty ±10 MeV. Note that T and vt parameters are anti¬ 
correlated. A surprisingly good agreement (see Fig. 1) of our model with the data on 
m-T-spectra serves as a strong support of the hypothesis of statistical nature of J/ip and 
ip' production [[TIJ and their kinetic freeze-out occurring directly at hadronization [Q. 

The dependence of the J/ip and ip' transverse mass spectrum on the centrality (quan¬ 
tified by the neutral transverse energy E T ) of Pb±Pb collisions at 158 A-GeV was also 
measured by NA50 Collaboration ||. An increase of (p?) and (p^) from peripheral colli¬ 
sions to the most central collisions can be explained, within our approach, by an increase 
of the model parameter vt with E T . Note that the increase of mean flow velocity and 
consequently an increase of (px) and (p\) with a centrality of the collision as well as an in¬ 
crease of (pt) and (p\) with particle mass m, are characteristic features of hydrodynamics. 
In contrast to J/ip and ip' mesons the my-spectra of the Drell-Yan pairs (dileptons with 
invariant mass M > 4.2 GeV/c 2 ) do not show this type of hydrodynamical behaviour. 
The values of (pt) and (p\) 0 for the Drell-Yan pairs are smaller than those for J/ip 
and ip' and do not change significantly with E 


The kinetic freeze-out parameters of pions were extracted from the analysis of single- 
and two-pion spectra measured for central Pb±Pb collisions at 158 A-GeV by the NA49 
Collaboration 
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21|1 . The results are: Tf = 120 MeV and vt — 0.55, i.e. they 


are very different than those obtained here from the analysis of heavy hadron spectra. 
In fact we checked that the parameters obtained for pions lead to the rrir spectra of 
heavy hadrons which strongly disagree with the data. A decrease of temperature and an 
increase of transverse flow velocity with time is a general property of expanding systems. 
The different kinetic freeze-out times of heavy hadrons and pions allow us to follow the 
expansion history of the hadron gas phase created in Pb±Pb collisions at 158 A-GeV. 
The freeze-out points of heavy hadrons and pions extracted from the data are plotted in 
Fig. 2 defining the path of the expanding hadron system in the T — TJt plane. 


In conclusion, the m^-spectra of D, J/ip and ip' produced in Pb±Pb collisions at 
158 A-GeV are analysed within the hydrodynamical model of the QGP expansion and 
hadronization. The spectra are in agreement with the hypothesis of kinetic freeze-out of 
these heavy hadrons occurring directly after the transition from the quark gluon plasma 
to the hadron gas. A mean collective transverse flow of hadronizing matter of vt — 0.2 is 
extracted from the fit to the spectra using temperature Th =170 MeV fixed by the analysis 
of hadron multiplicities [J3|. This result together with a previously obtained parameters 
of pion freeze-out (Tf = 120 MeV and vt — 0.55) allow for the first time to establish the 
history of the expanding hadron matter created in nuclear collisions. 

In the RHIC energy range the temperature parameter is approximately the same 
T = Th = 170 MeV ||, whereas the transverse hydrodynamic flow at T = Th is expected 
to be stronger. The model predictions for the mr-spectra of D and charmonia at the 
RHIC energies will be presented elsewhere. 
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Figure 1: The transverse mass spectra of (triangles down) and fl + (triangles up) [[HJ 
as well as J/-0 (circles) and ip' (squares) |J produced in Pb+Pb collisions at 158 A-GeV. 
The solid lines indicate a prediction of model (j^) assuming kinetic freeze-out of heavy 
hadrons directly after hadronization of expanding quark gluon plasma. The freeze-out 
parameters are: T = 170 MeV and vt = 0.194. 
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Pb + Pb at 158 A GeV 



Figure 2: The expansion history of strongly interacting matter created in Pb+Pb colli¬ 
sions at 158 A-GeV. The points indicate the temperature T and the mean transverse flow 
velocity Vt of matter at the time of J/rj; and i)}' freeze-out (upper point) and at the 
time of pion kinetic freeze-out (lower point). 



